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We consider a model where neutrinos interact with weakly interacting massive particles (WIMPs) 
as a possible explanation for the production of NuTeV's anomalous dimuon events. This model 
naturally produces the features reported by NuTeV, specifically the muon energy asymmetry and 
other kinematics. When combined with the NuTeV event selection criteria, this model predicts 
an expectation of zero events as observed in their other search channels. For the best agreement 
with data, the model requires a lower WIMP mass than conventionally expected, but one that is of 
current interest due to DAMA/LIBRA results, at around 10 GeV/c^. This model has a very high 
cross section for neutrino- WIMP interactions. We make suggestions for future searches that would 
confirm or refute this model. 

PACS numbers: 12.60.Jv, 13.15.+g, 14.80.-j, 95.35+d 



In the NuTeV experiment, a high energy {E « 
100 GeV) neutrino beam was directed at an instrumented 
iron target-calorimeter [H The primary goal of the 
experiment was to precisely measure the Weinberg angle 
sin^^u, J3j]. One of several other experimental goals in- 
cluded searching for neutral heavy leptons (NHLs). To 
enhance the sensitivity for NHL searches, a dedicated de- 
cay volume (decay channel) was introduced between the 
primary target (where neutrinos were produced) and the 
target-calorimeter (where the neutrinos interact). 

NuTeV collided 800 GeV protons with a beryllium ox- 
ide target. Secondary particles from the collision were 
steered by the magnets of the 'Sign Selected Quadrupole 
Train' (SSQT) to a decay region. The pious and kaons 
decayed to produce muons and neutrinos. The sign se- 
lection of the SSQT allowed for a very pure beam of z/^ 
(neutrino mode) or F^. The pion and kaon decay region 
was about 350 m long and was followed by a 900 m long 
shield to absorb everything but the i/'s. The decay chan- ' 
nel and the main detector followed. The decay channel 
consisted of a scintillator veto at the upstream end and 
drift chambers for tracking. Between the drift chambers 
it was filled with vessels of helium to minimize neutrino 
interactions. Details of this decay channel and its perfor- 
mance have been published [J-|6|. 

One analysis of the NuTeV decay channel data 
searched for final states of /i/i, fie, and /i7r. The ex- 
pected background of all three channels is 0.069 ± 0.010 , 
events, however there were three /i/i events in the data. 
These three events do not match characteristics of NHL's 
or any standard model sources. The most viable model 
proposed to date [7, 8] presumes that a neutral massive 
quasi-stable particle was created at the neutrino target ^ 
and decayed in the decay channel. We present an alter- 
native model here. 

Our model is based on the common speculation that 
a supersymmetric (SUSY) 0-[l2] ncutralino is a weakly 
interacting massive particle (WIMP) that accounts for ( 
dark matter. Dark matter is hypothesized to explain 




FIG. 1: The Feynman diagram for the stand-alone MC. 



discrepancies in galactic rotation curves and other astro- 
nomical anomalies Our hypothesis is that the 
neutrinos in the NuTeV beam interact with neutralinos 
from the dark matter halo present in the decay chan- 
nel. Charged current interactions excite the ncutralino 
(x°) to a chargino (x"*") state and produce a primary 
muon. When the chargino decays, it can produce a sec- 
ond muon as shown in Fig[TJ An initial study with Mad- 
graph 



171 revealed that there are hundreds of channels 



where , 



i + X ■ Madgraph currently 



assumes the initial particles are massless which is not a 
valid assumtion for this interaction. 

We developed a simple stand-alone Monte Carlo (MC) 
for the process shown in Fig. [T] Our analysis focuses on 
neutrino mode running because all three events occurred 
while NuTeV took data in that mode. Our model as- 
sumes an infinitely short lifetime so as to not spoil the 
vertexing. The expected lifetime of is very short and 
the relativistic effects should be minimal so this assump- 
tion seems reasonable. In a given execution of the MC, 
a specific mass for x° (M) and the mass difference be- 
tween it and (AM) were used. The x° was assumed 
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FIG. 2: Input distributions for the MC. The a) x and b) y 
locations, c) the neutrino energy, and d) the polar angle of 
the muon in the center of mass frame of the collision. 

30 

to be at rest. While this is not generally true, the initial 
momentum is expected to be negligible. The interaction 
location in z was taken to be uniformly distributed, while 
the neutrino energy and the x and y interaction profile 

5 were chosen to agree with the NuTeV larger data set of 35 
charged current events. These inputs are shown in Fig. 
[2I The neutrinos are assumed to have no transverse mo- 
mentum. In addition. Fig. [2Jl shows the polar angle of 
the primary muon in the center of mass frame of the col- 

10 lision. This angle was extracted from Madgraph as an 40 
input to the MC. 

Monte Carlo samples were generated at a variety of 



A few similar were tried using different variables. 
Specifically, the muon energy asymmetry was included in eo 

45 some definitions of x^, and some x^'s did not include the 
quantities tot or z. The set in Eq. 1 seems to include all 
variables with strong sensitivity and no varriables that 
are strongly correlated to others. By quantifying the 
agreement in this way, we were able to generate MC sam- 65 

50 pies and scan a large range of values for the masses of x*^ 
(M) and AM. We found that the NuTeV data gave a 
minimum x^ of 3.47 with 4 degrees of freedom at M — 
9 GeV/c^ with AM = 6 GeV/c^as you can see in Fig.H 
The range of scanned x° mass was from 0.5 GeV/c^ to 70 

55 1500 GeV/c^; while, AAf was scanned from 1 GeV/c^ to 
30 GeV/c2. 

The NuTeV events show reasonable agreement with 
the variables used in the x^ as shown in Fig. The Ef^a 



values of M and AM, and these samples were subjected 
to the same selection criteria used by the NuTeV anal- 
ysis. The selection criteria required exactly two muons 
with each carrying at least 2.2 GeV of energy, the sum 
of the energy of the two must exceed 12 GeV (to match 
the requirement in the other channels that the e or tt 
exceed 10 GeV), the angle of each muon was required 
to be less than 100 mr with respect to the beam direc- 
tion (ie the z-axis) , the muons needed to originate within 
the fiducial volume of the decay channel, be separated 
from the location of drift chambers, and each muon must 
reach the target calorimeter inside a transverse ±1.27 m 
fiducial area centered on NuTeV's coordinate system. Fi- 
nally, there was a transverse mass cut ttlt > 2.2 GeV/c 
where itit = \Pt \ + \J Pt + rnl where Pt is the dimuon 
momentum transverse to the beam direction and is 
the invariant mass of the visible particles. Reasonable 
agreement was observed between the NuTeV data and 
the simulation after event selection criteria were applied. 
This agreement will be discussed and shown shortly. 

In order to quantify the agreement between MC sam- 
ples and data, we defined a x^-like metric which used 
four input distributions: £"^1, £"^2, "^T, and z of the in- 
teraction. Muon number one is identified by its nega- 
tive charge when available and the higher energy is used 
when its charge is not available. For each quantity, the 
X^ compared the mean of the MC events passing selec- 
tion criteria to the mean of the three NuTeV events and 
treated the RMS of the NuTeV events as an uncertainty. 
The is defined as 



(1) 



shape is due to two categories of events being accepted by 
the selection criteria. In the rest frame of the W+, only 
a /i"*" going approximately parallel or anti-parallel to the 
rest/lab boost direction will be accepted by the selec- 
tion criteria. Those muons with a substantial transverse 
component have a negligible chance of being accepted. It 
is possible that this structure may allow additional sen- 
sitivity compared to using just the mean; however, we 
did not explore this possibility. The NuTeV events also 
show agreement with other kinematic variables in Fig [5l 
Specifically, NuTeV events sit firmly in the bulk of the 
distribution for the transverse radius of the interaction, 
the muon asymmetry (A= ^^^i,^^^^ ), the total px of the 
dimuon pair and invariant mass of the dimuons. 

The acceptance of the selection criteria applied to this 
model is very low. At the x^ minimum, the acceptance 
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FIG. 3: In the region of mass (M) and AM space with 
the lowest i contours of equal Sive shown. 
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FIG. 4: The input distributions for the its minimum 

(M = 9 GeV/c^ with a AAf = 6 GeV/c^) and at a lower AM 
value (M=2 GeV/c^ with a AA'/ = 3 GeV/c^). The arrows 
show the values of the 3 NuTeV events. The four distributions 
are a) z position, b) E^i of the negative muon, b) rriT and d) 



is 2.1 X 10~^, growing by about an order of magnitude 
as the mass of the x° decreases to 2 GeV/c^ and getting 
smaller as the mass increases. This behavior is quali- 35 
tatively expected: the larger the mass, the larger the 
possible opening angle between the muons, and the less 
likely both muons will hit NuTeV's main detector. This 
effect is also responsible for the steep shape of z for the 
accepted MC events, as shown in Fig |4^. The low accep- 40 
tance also means that, if this is an accurate model, other 
signatures should be observable. These will be discussed 
later. 

One feature of this analysis, which is unexpected in 
WIMP interactions, is the cross section. One would ex- « 
pect a cross section within a few orders of magnitude of 
lO"**^ cm^. However, when one considers the relatively 
small (10^^) number of neutrinos, the number of WIMP 
targets (10^) and the poor acceptance, in order to gen- 
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FIG. 5: Four other kinematic distributions where the x is 
at its minimum (M^o=9 GeV/c^ with a AM of 6 GeV/c^) 
and at a lower AM value (M = 2 GeV/c^ with a AM = 3 
GeV/c^). The arrows show the values of the three NuTeV 
events. The four distributions are a) the transverse radius of 
the interaction, b) the muon energy asymmetry, c) the total 
Pt of the dimuon pair and d) invariant mass of the dimuons. 



erate the three events in NuTeV, we would arrive at a 
cross section of around 10~^° cm^. While this is one of 
the least agreeable features of our model, it does make 
independent confirmation or rejection very feasible. 

In addition to the cross section, the fact that low- 
mass charginos and neutralinos were not observed by the 
LEP experiments in their examination of the Z width 
raises questions with our model. Some authors (for ex- 
ample [19|) have found regions of SUSY space that may 
still allow a neutralino with a mass like ours. To this 
point, our analysis has not made use of any SUSY spe- 
cific characteristics. So, to minimize these concerns, we 
will generalize our model. Our charged state is not neces- 
sarily a chargino, and our neutral state is not necessarily 
a neutralino, and, possibly, the mediator of the charged 
current reaction is not the standard model W. To mini- 
mize confusion throughout the rest of this paper, we will 
call the neutral state C°, the charged state C*, and the 
mediator D*. The C° and cannot interact with the 
Z boson with standard couplings. 

Dark matter candidates in the mass range that 
NuTeV's data prefer in our model have received a lot of 
attention lately (for example [20.425} and many of their 
citations) in response to the DAMA/LIBRA [2g| results. 
While they have interpreted the strong annual modula- 
tion of their events as dark matter interactions, other 
experiments [l^] have excluded their results. A feature 
of crystals, called channeling [28], has led some [29l. Iso!] 
to reevaluate the results, finding a small mass region be- 
tween 7-11 GeV/c2 or 3-8 GeV/c^ (2^ for the WIMP. 

A final concern is charged particle production (C'^C") 
at LEP. Actual chargino production at LEP may occur 
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through s-channel production via the Z or 7 and through 
t-channel sneutrino exchange. However, these modes de- 
structively interfere 3l| . In our model, the is not 
allowed to couple to the Z. It is reasonable to suggest eo 
that our C°, or another as yet unhypothesized particle, 
could participate in a t-channel interaction and the de- 
structive interference with s-channel 7 production may 
be more complete, producing a low enough cross section 
(^0.3 pb) to escape detection. Pair-production of our 65 

may also have escaped detection at LEP because re- 
cent searches (for example 32|, |33| ) for charginos are only 
sensitive above a AM of at least 3.0 GeV/c^. A reason- 
able portion of our favored space is available with AAf 
below 3.0 GeV/c^. For the simulated sample at M = 9/0 
GeV/c^ and AM = 3 GeV/c^, the distribution of E^a 
values are always below 3.5 GeV, in stark disagreement 
with two of the NuTeV events. The simulation at M — 2 
GeV/c^ and AA/ = 3 GeV/c^ has a good of 7.7 and a 
reasonable visual agreement. The distributions from this 75 
sample are shown in Figs. |4] and O 

Because multiple options may have allowed our model 
to have gone unnoticed until now, and, given the cur- 
rent interest in dark matter of order 5 GeV/c^, we have 
examined our model for ways in which the NuTeV col- so 
laboration may be able to confirm or refute this hypothe- 
sis. Instead of simulating backgrounds, we suggest event 
topologies for which the acceptance was relatively large 
and where we expect the backgrounds to be reasonable. 
One might think that NuTeV's examination 3J| of low as 
hadronic energy dimuons in their main detector should 
have seen events caused by neutrino-neutralino scatter- 
ing, however a few factors compete with better detector 
coverage. In that analysis, the lower energy muon was 
required to have an energy greater than 5 GeV and a « 
positive z momentum (Pz). In addition, NuTeV's main 
detector has a smaller fiducial volume. These lead us to 
expect fewer than 10 events due to our proposed model. 
If the energy requirement on this channel were lowered 
to 2.2 GeV, about 3 times as many signal events would 95 
be expected. 

There are a number of ways NuTeV could loosen their 
requirements and see more events in the decay channel of 
their detector. One option would be to require two tracks 
to be reconstructed, but allow one track to miss the de-100 
tector. The track that misses the detector would still 
be required to hit two drift chambers and reconstruct a 
good vertex. With these requirements, our model would 
expect between 20 (at low mass) and 45 (at high mass) 
events. It turns out NuTeV did not use the region lessios 
than 19 cm from the face of the main detector for the 
previous analysis. Presumably, this was in order to have 
multiple drift chamber hits from which to form tracks 
before the main detector (primarily for electrons and pi- 
ous). A second option, since we are interested in muonsno 
that can be found in the main detector, is to search, 
only the region from 19 cm to 2 cm from the face. For 



this search our model expects about 10 events with both 
muons reconstructed in the main detector. If we also re- 
move the requirements that both muons have an angle 
of less than 1 nirad we would expect around 25 events. 
A third option would be to remove the requirements on 
Tot, the angle and the energy of the second track, but still 
require both tracks to reach the main detector. NuTeV 
almost certainly required the muons energy above 2.2 
GeV to reduce standard model backgrounds when they 
were looking for a couple of events. However, with this 
selection we would expect at least 17 neutrino-neutralino 
scattering events. 

A clear signature of these events would be a pair of 
muons where one is a backward going muon in the lab 
frame. Backward going (P^ < 0) muons in the decay 
channel need to be restricted to regions where each muon 
hits at least two drift chambers, which produces a much 
smaller fiducial length and the acceptance is not much 
larger than their previous result. However, in the main 
detector, backwards going muons are an interesting op- 
tion. The backward going muons have very low energy, so 
good acceptance requires lowering the energy threshold 
to 1.1 GeV, which is still enough for significant pene- 
tration. If we look in the main detector for one muon 
with an energy greater than 9 GeV moving in the for- 
ward direction, one muon with an energy greater than 
1.1 GeV and a z component of momentum less than -0.02 
GeV/c, and a dimuon energy greater than 12 GeV, then 
we only expect one neutrino-neutralino scattering event 
at M=2 GeV/c^. But, by M = 6 GeV/c^ the expecta- 
tion has risen to 140 events and it stays at this level or 
above all the way to M = 75 GeV/c^. If backgrounds re- 
quire stronger constraints, one could insist that the back- 
ward muon penetrate at least four drift chambers (0.8 
m of iron) and the forward muon have it's momentum 
and charge analyzed in the toroid. Then, we expect no 
events at very low mass (M = 2 GeV/c^) but for M — 10 
GeV/c^ and above we expect 80 neutrino-neutralino scat- 
tering events. 

In addition to NuTeV, a number of past, current, and 
future experiments should be able to search for similar 
topologies. At our lowest , the neutrino energy thresh- 
old is 10 GeV, so MINOS or MINERz/A would likely only 
be able to study this if they ran the NUMI beam in high 
energy mode. CHORUS and NOMAD both have high 
enough energy, should have enough inte nsity , and have 
looked at dimuon production before [35l |36|. NOMAD 
has studied backward going particles [37|. Given these 
experiences, they should be able to look for clean two 
muon events with one muon going backwards and per- 
haps quickly rule out this explanation for the NuTeV ex- 
periment's dimuon excess. In the future, if the NuSoNG 
[38| experiment were brought to life, they would expect 
60 events with the same topology. By instrumenting the 
sides of their helium region, they ought to be able to 
quickly confirm or refute this explanation. 
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One might also ask why the current cohider detectors 
have not noticed signs of new physics due to this model. 
It is possible they already have. In Refs. [H, li^, the 
CDF collaboration reports observing a significant num- 
ber of relatively low momentum muons which appear to 
originate very far from the interaction point. If our 
has a lifetime somewhat larger than a B meson's, then 
pair production of would appear this way. Alterna- ' 
tively, b and c quarks produced in interactions, would 
decay producing neutrinos which would interact with 
dark matter in the detector volume. However, the DO 
collaboration searched and was unable to confirm this < 
signal[41]. In addition to pair production of the charged 
state, hadron colliders could produce a mixed state with 
C^C°. That assumes, however that the D^ interacts 
with quarks. Pair production of the neutral state as- 
sumes some Z-like object which interacts with the neutral 
state and quarks. If these final states can be produced, 
we would need a value for the coupling to quarks in order 
to compare to the Z -I- jet and /i + MET topologies that 
will conceal them. ' 

Our hypothesis has a number of attractive features in 
that it provides an explanation for the NuTeV dimuon ex- 
cess, and the kinematics, especially the asymmetry, are 
in reasonable agreement with NuTeV's events. Our hy- f 
pothesis is also in agreement with the feature that all 
of NuTeV's events were in neutrino mode, because the 
integrated neutrino beam had about 3 times as much in- 
tensity as the anti-neutrino. It is a natural result of the 
selection criteria in this model that there will be no fie 
and no fiir events because the flavor of the high energy 
lepton is determined by the incident neutrino flavor, and 
the particle due to the decay is always below the e or 
TT threshold of the NuTeV analysis. In addition, this hy- - 
pothesis has the interesting feature that it indicates a low 
mass dark matter candidate. 

We have considered a simple model to explain NuTeV's 
anomalous dimuon events. In this model, the neutrinos c 
undergo charged current scattering off dark matter par- 
ticles and excite them to a more massive charged state. 
When the massive charged state decays back into the 
dark matter particle, the net visible signature is two 
muons. This model nicely describes the kinematics and 
other features of the NuTeV events, and those kinematics 
are used to extract a value for the mass of the dark matter 
particle and the mass of its charged partner. The most 
unsettling features of this model are the cross section," 
which is many orders of magnitude above expectations, 
and the unobserved low mass charged particle, C^. We 
have provided a number of suggestions for further studies 
which could confirm or reject this explanation. 
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